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Abstract 

We consider the two point correlation function of scalar current in QCD 
sum rules approach to estimate the overlap amplitude of oq meson. We then 
employ QCD sum rules to calculate the coupling constants ga;ao7 aiid gpao7 
by studying the three point 00^7- and ao/97-correlation functions. 
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I. INTRODUCTION 

The low-mass scalar mesons have fundamental importance in understanding the theory 
and phenomenology of low energy QCD. From the experimental point of view, isoscalar 
/o(980) and isovector ao(980) are well established, but the nature and the quark substructure 
of these scalar mesons, whether they are conventional qq states |I[], KK molecules [0], or 
multiquark exotic q^lf states ||^, have been a subject of controversy. On the other hand, 
since they are relevant hadronic degrees of freedom, besides the questions of their nature, 
the roles of scalar mesons in the hadronic processes must be studied. 

The radiative decay processes of the type V^ — » P^P^'y where V and P belong to the 
lowest multiplets of vector (V) and pseudoscalar (P) mesons have become a subject of re- 
newed interest because they offer the possibility of investigating the new physics features 
governing meson physics in the low energy region. Although these rare decays have small 
branching ratios due to absence of bremsstrahlung radiation, their study offers the pos- 
sibility of testing the theoretical ideas about the interesting mechanisms of these decays, 
as well as shedding light on the structure of intermediate states involved in these decays. 
Particularly interesting are the exchange mechanisms of scalar resonances contributing to 
these decays. The radiative decays p° -^ tt^i]'-/ and uj -^ Tf^rj'-f were studied using a low 
energy effective Lagrangian approach with gauged Wess-Zumino terms [Q, and later by us- 
ing standard Lagrangians obeying SU(3)-symmetry [^. In both of these calculations scalar 
meson intermediate state contributions were neglected and the contributions of intermediate 
vector mesons were taken into account. However, it is of interest to study the contribution 
of ao-intermediate state to these decays as well, and for that a knowledge of 00^7- and 
aop°7-vertexes are needed. 

In this work, we estimate the coupling constant gaop7 and gaooj-y by employing QCD sum 
rules which provide an efficient method to study hadronic properties and which have been 
employed to study hadronic observables such as decay constants and form factors in terms 
of nonperturbative contributions proportional to the quark and gluon condensates 
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II. ANALYSIS AND RESULTS 

The QCD sum rules approach [^^ is a model independent method to study the prop- 
erties of hadrons through correlation functions of appropriately chosen currents. We choose 
the interpolating currents for uj and p mesons as j|^ = |(m7M + d-^d) and j^ = \{u'^u — d'~^d) 
respectively; and for oq meson as jao = \iuu — dd) [§|,0], and we work in the SU(2) flavour 
limit m„ = m^ = niq. In the sum rule, the overlap amplitude of a^ meson Xa^ =< O|jaoko > 
is needed. In a previous work [^] we studied the scalar-isoscalar a meson by considering the 
two-point scalar current correlation function. Since perturbative and QCD-vacuum conden- 
sate contributions to scalar current correlation functions cannot distinguish between isoscalar 
and isovector channels, we follow here the same method and we study the scalar-isovector 
Oq meson by considering the two-point current correlation function 

n(p2) = ^|dW^-<o|ra,(x)i(o)}|o> . (1) 

The two-loop expression for the scalar current correlation function 11 (p^) in perturbative 
QCD was calculated |T^, and for light quark systems in the limit niq = it is given by the 
expression 
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QCD-vacuum condensate contributions to the scalar current correlation function 11 (p^) were 



obtained by the operator product method ||Tl| in the same limit mg=0 as 



Let us note that the term < mqqq > is independent of quark mass since it is given as 
—f^ml/A through Gell-Mann-Oakes-Renner relation 0. 

The correlation function n(p^) satisfies the standard subtracted dispersion relation [^ 

n,e.t(/)=p'/ -, ^p(s) + n(o) (4) 

JO s[s — p'') 



where the spectral density function is given as p{s) = {l/7i)l'mll{s). The spectral density 
contains a single sharp pole irXaoS^s — rnl^) corresponding to the coupling of a^ meson to 
the scalar current. The continuum contribution of the higher states to the spectral density 
is estimated as p = ph{s)9{s — sq) where sq denote the continuum threshold and ph is given 
by the expression Ph{s) = {l/n)ImIlopE{s) with IIope{s) is obtained from Eq. (2) and Eq. 
(3) as IIope{s) = Ilpert{s) + Ucondis). After performing the Borel transformation we obtain 
the QCD sum rule for the overlap amplitude Xao as 
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+ 2M^ < -«^^ > +Y^^2 < »sG > -^^ < a.m > . (5) 

In the numerical evaluation of Eq. (5) we use < rriqqq >= (—0.82 ± 0.10) x lO^'' GeV*, 
< a,G2 >= (0.038±0.011) GeV\ < as{qq? >= -0.18 x lO^^ GeV^ ^^. The threshold 
is choosen below a possible ao(1450) pole and it is varied between so=1.6-1.7 GeV^. Since 
the Borel parameter has no physical meaning, we look for a range of its values where the sum 
rule is almost independent of M^, we choose the interval of values of Borel parameter M^ as 
1.2-2.0 GeV^. The overlap amplitude Xa^ as a function of M^ in this interval for different 
values of sq is shown in Fig. 1 from which by choosing the middle value M^=1.6 GeV^ in it 
is interval of variation, we obtain the overlap amplitude as A^o = 0.21 ± 0.05 GeV^ where 
we include the uncertainty due to the variation of the continuum threshold and the Borel 
parameter M^ as well as the uncertainty due to errors attached to the estimated values of 
condensates as quoted above. 

In order to derive the QCD sum rule for the coupling constants gaou-y and gaop'y, we 
consider the three point correlation function 

t,Ap,p) = J d'xd'ye'p'-ye-'"-^ < o|r{j;(o)j:(x)j„„(i/)}|o > , (6) 

where j^ = {cyU'-y^u + e^d'y^d) is the electromagnetic current with e„ and Cd being the quark 
charges, and j^ is the interpolating current for u or p° meson. 



In order to obtain the phenomenological part of the sum rule, we consider the double 
dispersion relation for the vertex function T^^, 

where the possible subtraction terms are neglected since they will not make any contribution 
after double Borel transform. For low values of Si and S2i the spectral function p^,y(si,S2) 
contains a term proportional to double (5-function 5{si — m\r)5{s2 — iti'^^^), corresponding to 
the transition ao -^ V'~f where V denotes uj or p^ meson. We therefore saturate the dispersion 
relation satisfied by the vertex function T^j, by these lowest lying meson states in the vector 
and the scalar channels, and this way we obtain for the physical part 

rj. f ,, _ < 0\j^\V X V{p)\j;\ao{p') X aoljaolO > , 

where the contributions from the higher states and the continuum is denoted by dots. In 
this expression the overlap amplitude \a^^ =< ao|jao|0 > of ao-meson has been determined 
in previous sections. The overlap amplitude Ay of vector meson is defined as < 0|j^|V >= 
XyUi, where Ui, is the polarization vector of the vector meson u or p°. The matrix element 
of the electromagnetic current is given as 

<V{p)\Jl\ao{p')>= -i 9aoV^K{q^){p-qUf,-u-qPf,) , (9) 

mv 

where q = p — p' and K{q^) is a form factor with K(0) = 1. This expression defines the 
coupling constant through the effective Lagrangian 

C = —ga^v^d'^V^idaAfs - dfsA^)ao (10) 

mv 

describing the aoV"7-vertex. 

The theoretical part of the sum rule is obtained by calculating the perturbative contri- 
bution and the power corrections from operators of different dimensions to the three point 
correlation function T^j,. For the perturbative contribution we consider the lowest order 
bare-loop diagrams shown in Fig. 2(a). Furthermore, we consider the power corrections 



from the operators of different dimensions that are proportional to vacuum condensates 
<qq>, <qa ■ Gq > and < (qq)^ >■ We do not consider the gluon condensate contribution 
proportional to < G^ > since it is estimated to be negligible for light quark systems. We 
perform the calculations of the power corrections in the fixed point gauge [0. We work 
in the limit rriq = 0, and in this limit perturbative bare-loop diagram does not make any 
contribution. Moreover, in this limit only operators of dimensions d=3 and d=5 make con- 
tributions that are proportional to < gg > and < qa ■ Gq >, respectively. The relevant 
Feynman diagrams for power corrections are shown Fig. 2(b) and (c). If we consider the 
gauge invariant structure {Pfj,qu — P ■ qOfiu), we obtain the power corrections of dimensions 
d=3 and d=5 as 

3 11 - 

C3 = i- — —^{eu <uu> +ed< dd>) (11) 

and 
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After performing double Borel transform with respect to the variables Q^ = —p^ and 
Q' = —p' , and by considering the gauge-invariant structure (p^g^ — p ■ qg^iu) for the phe- 
nomenological part as well, we obtain the sum rule for the coupling constant gaoy7 

3my <) :4 /3 9 2 1 1 2 1 



g^,y, = -e, <uu> J^e^e^ (^- - -m,— - ^m,^) (13) 

where e^ = {cu + e^) for p° meson and Cg = (e„ — e^) for u meson, and we use the relations 
< qa ■ Gq >= m^ < qq > and < uu >=< dd >. For the numerical evaluation of the 
sum rule we use the values mg = (0.8 ± 0.02) GeV^, < uu >= (-0.014 ± 0.002) GeV^ 



PJT^, and nip = 0.770 GeV, m^ = 0.782 GeV. For the overlap amplitude A^q we use 
the value Xao = (0.21 ± 0.05) GeV^ that we have estimated previously. We determine 
the overlap amplitude Ay for uj and p° meson from the measured leptonic decay widths 
T(y -^ e^e^) |]r^, thus we use their experimental values Ap = (0.117 ± 0.003) GeV'^ 



and A^^ = (0.108 ± 0.002) GeV"^. In order to analyze the dependence of gagy^ on Borel 



parameters M^ and M' , we study the independent variations of M ^ and M' in the interval 
0.6 GeV^ < M'^,M' < 1.4 GeV^ as these hmits determine the allowed interval for the 



vector channel [|I^. The variation of the coupling constant gajao7 as a function of Borel 
parameter M^ for different values of M' is shown in Fig. 3, examination of which indicates 
that it is quite stable with these reasonable variations of M^ and M' . We choose the 
middle value M^ = 1 GeV'^ for the Borel parameter in its interval of variation and we 
obtain the coupling constant gaouj^ as gaocu-y = (0.75 ± 0.20). We indicate the error arising 
from the numerical analysis of the sum rule as well as from the uncertainties in the estimated 
values of the vacuum condensates. In Fig. 4 we present the variation of the coupling constant 
gaopj as a function of the Borel parameter M^ for different values of M' . Following a similar 
analysis as in the case oi ga^u-y, we obtain the coupling constant gaop7 as gaop7 = (2.00±0.50). 
The values for the coupling constants gaocj7 and gaop7 that we obtain are in agreement with 
the expected SU(3) ratio gaQp-y:gaQ^^=3:l. 

In our analysis, we use for the overlap amplitudes A^; and Xp the values that we obtain 
from the experimental electronic decay widths of u; and p° mesons. On the other hand, it may 
be argued that in a QCD sum rules calculation it is more appropriate to use the values of the 
overlap amplitudes that are also determined within the framework of QCD sum rules method. 



Electromagnetic decays of vector mesons using QCD sum rules were studied [0, and in their 
analysis the authors used the values of the overlap amplitudes A^^ = (0.16 ± 0.01) GeV^ 
and Xp = (0.17 ±0.01) GeV^ that they also determined utilizing QCD sum rules. If we use 
instead these values of the overlap amplitudes in our calculation we then obtain the coupling 
constants gaow7 and gaop7 as gaow7 = 0.45 ±0.10 and gaop7 = 1.30 ±0.30 which are consistent 
with our above results. 

In the investigations of the role of oq meson in hadronic processes, the relevant coupling 
constants of oq meson are needed. In this work, we employed QCD sum rules approach to 
estimate the coupling constants gaoi^7 and gaoP7- We feel that the studies of the different cou- 
pling constants of oq meson should be continued. In particular QCD sum rules calculations 



should be improved by taking into account the high order corrections to the perturbative 
part of the three point correlation function and also to the two point correlation function 
employed in the estimation of the overlap amplitude of oq meson. 



REFERENCES 

[1] N. A. Tornqvist, M. Ross, Phys. Rev. Lett. 76 (1996) 1575. 

[2] J. Weinstein, N. Isgur, Phys. Rev. D 41 (1990) 2236. 

[3] R. L. Joffe, Phys. Rev. D 15 (1977) 267; D 17 (1978) 1444. 

[4] S. Fajfer and R. J. Oakes, Phys. Rev. D 42 (1990) 2392. 

[5] A. Bramon, A. Grau and G. Pancheri, Phys. Lett. B 283 (1992) 416. 

M. A. Shifman, A. L Vainstein, V. L Zakharov, Nucl. Phys.B 147 (1979) 385 and 448. 

[7] L. J. Reinders, H. R. Rubinstein, S. Yazaki, Phys. Rep. 127 (1985) 1. 



P. Colangelo, A. Khodjamirian, liep-ph/0010175 , to be pubhshed in the Boris loffe 



Festschrift, World Scientific, Singapore, 20001. 
[9] A. Gokalp, O. Yilmaz, Phys. Rev. D 64 (2001) 034012. 
[10] K. G. Chetyrkin, Phys. Lett. B 390 (1997) 309. 
[11] E. Bagan, J. L La Torre, P. Pascual, Z. Physics. C 32 (1986) 43. 
[12] B. L. loffe, A. V. Smilga, Nucl. Phys. B 232 (1984) 109. 
[13] A. V. Smilga, Yad. Fiz. 35 (1982) 473 [Sov. J. Nucl. Phys. 35 (1982) 271]. 

[14] V. M. Belyaev, B. L. loffe, Zh. Eksp. Teor. Fiz. 83 (1982) 876 [Sov. Phys. JETP 56 

(1982) 493]. 

[15] Review of Particle Properties, D. E. Groom et al., Eur. Phys. J. C15 (2000) 1; 
[16] V. L. Eletsky, B. L. loffe, Ya. L Kogan, Phys. Lett.B 122 (1983) 423. 
[17] Shi-lin Zhu, W-Y. P. Hwang, Ze-sen Yang, Phys. Lett. B 420 (1998) 8. 



Figure Captions: 

Figure 1: The overlap amplitude Aqj, as a function of Borel parameter M^. 

Figure 2: Feynman Diagrams for the 00^7- vertex: a- bare loop diagram, b- d=3 operator 
corrections, c- d=5 operator corrections. The dotted lines denote gluons. 

Figure 3: The coupling constant Qa^ui-y as a function of the Borel parameter M^ for different 
values of M' . 

Figure 4: The coupling constant Qaofrt as a function of the Borel parameter M^ for different 
values of M' . 
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Figure 4 
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